Theory of Multi-Band Superconductivity in Iron Pnictides 
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The precise nature of unconventional superconductivity in Iron Pnictides is presently a hotly 
debated issue. Here, using insights from normal state electronic structure and symmetry arguments, 
we show how an unconventional SC emerges from the bad metal "normal" state. Short-ranged, multi- 
band spin- and charge correlations generates nodeless SC in the active planar dxz,yz bands, and an 
inter-band proximity effect induces out-of-plane gap nodes in the passive d^^2_^2 band. While very 
good quantitative agreement with various key observations in the SC state and reconciliation with 
NMR and penetration depth data in the same picture are particularly attractive features of our 
proposal, clinching evidence would be an experimental confirmation of c-axis nodes in future work. 

PACS numbers: 74.70.-b, 74.20.-z, 71.30.-|-h, 74.20.Rp 



High Temperature Superconductivity (HTSC) in the 
recently discovered Iron pnictides (FePn) is the latest 
surprise among a host of others in d- and / band mate- 
rials [ij. HTSC arises upon doping (x) a spin-density- 
wave (SDW) metal: with x > Xc, SDW order melts, 
giving way to a HTSC. In the "normal" phase, with no 
SDW/SC order for T > Tsjjyy/sc{^)t various probes re- 
veal an incoherent metallic state, presumably not too 
far from a Mott insulator HTSC in FePn arises 

from an incoherent normal state, at least in the FeAs- 
based materials, where a quasi-linear-in-T resistivity is 
seen immediately above T^. this implies a correlated 
normal state. Local density approximation plus dy- 
namical mean-field (LDA-I-DMFT) calculations in the 
intermediate-to-strong correlation limit [^, indeed find 
very good quantitative agreement with a range of basic 
physical properties of FcPn, attesting to their sizably cor- 
related nature. 

//SC is an instability of this incoherent metal, one ex- 
pects pronounced deviations from weak-coupling BCS- 
like theories in the SC state, as in the cuprates: this 
is indeed borne out in experiments. Specifically, (i) the 
pair coherence length, S^ab is short (< 20 A), the up- 
per critical field, iJc2 is high with small superfluid 
density, (ii) the NMR relaxation rate, Tf ^(T) ~ 
and the penetration depth in the 1111-FePn, A(r) ~ T" 
with n = 2 — 2.5 7|, indicating unconventional SC (U- 
SC). However, the 122-FePn seem to show conventional 
s-wave-like behavior, though, even here, the situation is 
controversial 0]- Moreover, /zSR studies reveal Uemura 
scaling and Tc vs ep scaling characteristic of other, well- 
known, correlated SC 0]. (iii) Optical measurements [§| 
show a large-scale spectral weight transfer (SWT) across 
Tc over an energy scale > 0(1. 0) eV, while ARPES 
data show a sharpening of the low energy quasiparticle 
kink [13] below Tc- Taken together, as they must be, 
(i)-(iii) imply a strong coupling SC arising from an in- 
coherent non-FL metal. Given (ii), the existence of gap 
nodes is an open, controversial, issue [ll|. In particular, 
(i)-(iii) resemble observations in cuprates implying 



an U-SC closer to the Bose condensed, rather than the 
BCS limit [lo|. How can such a U-SC arise as an in- 
stability of an incoherent normal state [3, How may 
nodes appear in the gap function? Does SC pairing in- 
volve all, or a subset of the Fe-d orbitals? These issues 
are of great import for FePn, and call for a systematic 
resolution. Here, motivated by the above, we focus on 
the 1111-FePn, specifically on LaOi-xFsAsF^. Guided 
by group-theoretical analyses, and by the renormalized 
multi-band structure, we propose a specific, strong cor- 
relation based, interband pairing mechanism for U-SC as 
an instability of the incoherent metal found in earlier 
LDA-hDMFT studies SQ. 

LDA-hDMFT studies for the 1111-FePn find an inco- 
herent, pseudogapped "normal" state, implying blocked 
coherent one-electron propagation: only collective mo- 
tion of the spin and charge fluid is possible [l3|- This 
opens the door to coherent fwo-particle propagation, i.e, 
to two-particle instabilities, to relieve the finite entropy 
at lower T. At a; = 0, the electron- and hole-like Fermi 
surface (FS) sheets are nearly nested in LDA calculations, 
and DMFT will not alter this fact, since the self-energies 
are local; this favors the q = (tt, 0) SDW instability, also 
within LDA+DMFT [31 (in fact, no nesting is needed 
at strong coupling |3||), in agreement with neutron 
data 0, [3| . Doping weakens the SDW and facilitates 
the U-SC instability. 

Since coherent one-electron transfer between different 
orbitals is blocked in the incoherent metal, we propose, in 
analogy with the situation encountered in coupled Lut- 
tinger liquids (LL) [15], that small, residual, inter-site 
and inter-orbital interactions mediate two-particle insta- 
bilities in FePn at low T. In principle, various types of 
instabilities, competing with U-SC [l6|, may be gener- 



ated, as studied more extensively for cuprates 17[. Here, 
we focus on the U-SC alone. We further assume that pair 
formation primarily involves the active dxz,yz orbitals, 
and show that, due to the inter-orbital coupling U' , an 
interband (IB) proximity effect will induce SC on all FS 
sheets, as in Sr2RuOi 18| . 
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Excluding spin-triplet pairing, the general interaction 
in the cooper channel is 



Hpair — r) ^afcC^' ^ )'^l,fc,t'^I-fc,i 



Cb,-k',lCa,k','\ , 



a^h.k.k' 



where a,b = xz,yz and the scattering vertex is 
Vab{k,k',u!) = g'^Xabik - k' ,uj) with Xab{k - k',uj) be- 
ing the mier-orbital susceptibility. Such a term natu- 
rally arises at second order from a one-electron inter- 
band term, tabJ2<i,j>,ai^ia<j^jbo- + h.c), when the one- 
electron spectral function is incoherent, much like in cou- 
pled LL [l9|. This is an important feature of our work: 
normal state one-particle incoherence, arising from the 
Anderson orthogonality catastrophe 3, 4] due to strong 
multi-orbital correlations, favors two-particle coherence. 
The static, nearest- and next-nearest neighbor parts of 
Vabik,k') are v}l\k,k',0) ~ " <^(40 - 50) meV 

0(15 — 20) meV, close to the su- 



(2) 



and V;, _ jjT+j^ 
perexchange scale estimated from inelastic neutron scat- 
tering (INS) studies ^2d\. Notice that the effective inter- 
action contains coupled inter-orbital charge and spin fluc- 
tuations involving nearest (n.n) and next-nearest neigh- 
bor (n.n.n) Fe sites, and is explicitly k-dependent, favor- 
ing U-SC from the outset. Hpair can be decoupled a la 
Gorkov 2l| to yield 



ttMF _ 



a,b.k 



[Aab{k)clJ 



^^-k,b,i 



where Aabik) = ^Vab{k){c^k,b,l'^k,a,i;) is the SC gap 
function. Consistent with lattice symmetry, and in- 
cluding the n.n and n.n.n contributions, Vab{k,k') = 
J2i^lbVi{k)vi{k'), where the r]i{k) are irreducible 
reprersentations of D4^h poin t g roup. So we can ex- 
pand Aabik) = Ei^Lmik) iJ. Explicitly, Aabik) - 
Aiicx + Cy) + A2CxCy, with Cq, =cosfcQ,. The case A2 = 
(-629 representation of D^h) is favored in the numerical 
study ^22l|, however, A2 7^ is rigorously required [23| . 
With well-separated electron- and hole-like FS sheets as 
in LDA, no in-plane nodes exist in the SC gap for A2/A1 
in the chosen range. A A2 ^ is also favored by the ob- 
servation of appreciable geometric frustration in FePn, 
where J2/J1 — 0(0.7 — 1.0) have been deduced from ab- 
initio studies, and consistent with INS results [20j. 

The "normal" state is modelled by a five-band 
Hubbard model, iJ„ = Hn + Hi, treated ear- 
lier with LDA+DMFT [1 |3|, and reads Ha = 

■ Si^p, where a, /S — 



T,k,a^k,c.cl Ck,a,a while Hi 



r . In the U-SC phase, we have 



xy, yz, zx, x'^ — y^, 3z^ 
to solve H = Hn + -ffpiv within LDA-fDMFT. For- 
tunately, the intersite nature of Hpair adds simply a 
term bilinear in fermions to H . The LDA-t-DMFT equa- 
tions are now readily extendable to the U-SC regime: 



the Gab and S^h now have normal and anomalous com- 
ponents, and a closed set of DMFT equations, yield- 
ing both Gaa,Fab selfconsistently, is solved by extend- 
ing earlier LDA-I-DMFT [3] to include an explicit pair- 
field term. Including the pair-field, the propagators are 



i^+iikb+^lii^) ^ and 
where the * denotes 



Gaaik,Lu) = [UJ - Cka - Sa(w) - 

Fabik,u:)^Gaaik,u:) ^^^;;^^,^^^ , 

complex conjugation. Given intersite pairing, Fij^abi^^) 
falls off at least as (here, D is the lattice dimen- 

sion) the dynamical effects of these non-local pair fluc- 
tuations do not enter the DMFT self-energies. Obvi- 
ously, however, they do affect the GaaS- Our treat- 
ment is thus different from that of Garg et al. fj^, 
where the dynamical effects of local pair fluctuations 
must be kept in DMFT. Anomalous contributions like 

^S,abit) = -V^,G,,,,aait)F,,,^abi-t)G,,,,bbit)^ etC, tO 

the dynamical self-energy would thus enter in a cluster- 
DMFT approach, but drop out in D = cxd. It then suf- 
fices to use the DMFT with the modified matrix prop- 
agators above. Finally, since these equations couple all 
Fe-d-orbitals, the opening up of a SC gap in the dxz,yz 
bands could induce secondary gaps in the remaining d 
orbitals, in a way reminiscent of the inter-band proxim- 
ity effect in Sr2RuO/i [3| (see below). Further, given 
sizable U = 4.0 eV, U' = 2.6 eV for FePn 0, [|, large 
spectral changes within DMFT-like approaches should 
accompany the U-SC instability: these are seen in spec- 
tral probes [1, @| . Is such a scenario also consistent with 
other signatures of a strong couphng, U-SC? 

Aiming to shine light on these issues, we now describe 
our results. Upon convergence of the DMFT equations, 
the one-electron- and pair spectral functions can be read 
off and used for direct comparison with observables in 
the U-SC state. In Fig [H we show the changes induced 
by U-SC in the total one-eletron DOS. Clear sharpening 
of the low-energy kink at ~ 20.0 meV^SUHli already 
somewhat visible in the "normal" state [26[ , is seen across 
the U-SC instability. Strong normal state incoherence 
iTiliuj) in the equation for Gaa(k, w)) prevents opening 
up of a clean SC gap in the DOS, as indeed observed 
in PES. Remarkably, the position of this kink and the 
detailed PES lineshape up to a binding energy of —0.5 eV 
are both quantitatively reproduced within our theory. 

More microscopic insight is obtained by investigating 
the orbital resolved DOS. As seen from Fig[2l only the 
dxz,yz and d3z2_r2 DOS show the sharpening of the kink 
feature. Moreover, dominant SWT occurs from the dxy 
band to the dxz,yz,3z^-r^ bands, as clearly seen in Fig. [21 
The first implies an orbital selective coupling of the car- 
rier propagators to multi-orbital, overdamped and short- 
ranged, charge- and spin correlations, and fully agrees 
with indications from ARPES studies [2^, which show 
clear evidence thereof. Large SWT, over an energy scale 
0(2.0) eV, also accompanies the U-SC instability. This 
is also in qualitative accord with findings in optical stud- 
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FIG. 1: (Color online): Evolution of the total one-particle 
density-of-states (DOS) across the superconducting transition 
for doped {x = 0.1) LaOi-^FeAsF^, for U = 4.0 eV, U' = 

2.6 eV, and the effective interaction, Vat =™£^x(V^'^' , Vj^') ~ 
50 meV. The inset shows how the low-energy kink, already 
visible above Tc (dashed curve) , sharpens up below Tc (dotted 
curve), in very good agreement with ARPES data [lol. [25j. 



FIG. 2: (Color online): Orbital resolved evolution of the 
one-particle DOS across Tc, for the same parameter set as 
in Fig. [T] The orbital-selective gap structure, arising from 
multi-band correlations, is manifest. This is manifested in the 
orbital-selective sharpening of the low-energy kink, in quali- 
tative accord with ARPES [13, [111. 



ies fs^. Taken together, these findings imply a strong 
coupling SC, and fully accord with other observations in 
the SC state in FePn, as shown below. 

Additional fundamental features of interest are gleaned 
from observation of the sharpening of the kink in the 
d^z^^r^ DOS in the U-SC state above. This strongly sug- 
gests that an IB-proximity driven, out-of-plane SC gap 
component should enter on symmetry grounds. Exami- 
nation of the Slater-Koster fit to the LDA shows a sizable 
d^z'^-r'^ mixing 27|; this will induce a A'cosffc^c) 



component, reduced by dipolar reduction factors [l8|, in 
the SC gap function. Strong supportive evidence for this 
comes from extant dHvA studies on LaFePO This 
indeed shows significant warping of the dxy band along 
the c-axis in the non-magnetic state, caused precisely by 
the mixing with the d^z^_j.^ pocket states. While cor- 
relation effects (LDA+DMFT) result in band shifts and 
dynamical SWT, they will leave the LDA FS unchanged, 
since the self-energies are purely local: this observation 
should hold for the 1111-FePn in general. Such interband 
one-particle mixing between the xy and 3z^— bands in- 
duces an additional cos(fczc) component in Aa;,(k), giving 
possible line nodes in the SC gap function along k^. This 
may be theoretically checked by closely following earlier 
analysis for another multi-band U-SC, Sr2Ru04 

Generation of line nodes along c is not in 
conflict with extant ARPES results, as seen in 
{Sr / Ba)i^xKxFe2As2 [l^. Moreover, it can naturally 
explain the NMR and penetration depth data, being ca- 



pable of giving the observed power-law-in-T behavior for 
both the NMR Tf ^(T) and A„b(T) ff\ at low T « T^. 
In contrast to the s± gap with disorder scenario [2^ . ex- 
istence of nodes would imply universal power-law depen- 
dences in these quantities as a function of x, and along 
differing FePn members, as recently seen Q. 

We emphasize that whether nodes along exist or 
not is a delicate matter, and that a careful examina- 
tion of the details of the band structure should pro- 
vide invaluable clues. For example, dHvA study [28| on 
LaFePO indeed finds sizable c-axis corrugation of the 
dxy band, indicating sizable IB-proximity, and, interest- 
ingly, \ab{T) ~ T^-^ [2I], strongly suggesting Hue nodes 
in the SC gap. However, data on SmFeAsOi-xFx (3^ 
are consistent with a smooth (nodeless) angular variation, 
A(0) = Ao(l + ecos4(/)), of the in-plane gap, and Xab{T) 
follows the well-known e~^('^)/'^^ law. In the FeAs mate- 
rials, it is conceivable that the multiband coupling does 
not always induce out-of-plane nodes in the SC gap; this 
will be sensitive to the detailed topology of the renormal- 
ized FS. In cases where a non-exponential A(r) is seen, we 
suggest that dHvA results should show significant c-axis 
warping of the dxy-FS sheet: this should be investigated 
in more detail. However, the planar (1 -I- ecos4</)) varia- 
tion deduced for Sm-based FePn is not inconsistent with 
our in-plane form factor, as proposed above. Addition- 
ally, recently, INS work on iVi-doped BaFe2As2 clearly 
shows both planar and c-axis variation of A(k); inter- 
estingly, precisely lS.{kz) ~cos{kzc) was found there [lH . 
It would be very interesting to see if this is also true 
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of the llll-FePn. We then predict that this should be 
correlated with the c-axis warping of planar d^y band in 
dHvA work. This would constitute a non-trivial check of 
our proposal. 

From LDA-t-DMFT, we estimate the gap magnitude 
to be Aab =i 4.5 meV, close to that found in PES [H]. 
Using the carrier Fermi velocity from LDA (note that 
this is not changed by LDA+DMFT), vp ^ 0.7 eV.A, 
the relation ^ ~ vf^/t^A gives ^ ~ 60 — 80A, close to 
observed values. The upper critical field is now estimable 
as Hc2 = 0o/27r^^ ~ 50 - 70 T, (0o is the usual flux 
quantum) again consistent with experiment [5|. These 
are classic signatures of a strong couphng SC [3j|, and 
require an instability of a strongly correlated incoherent 
"normal" state, as proposed here. 



family. Our proposal calls for careful experimental search 
for the out-of-plane nodes in the 1111 Iron Pnictides. 



Our proposal is distinct from other, extant ones [11 



El,!!^. The s± proposal [lil. [sot is derived from 
a weak-coupling instability of an itinerant FL. Wu et 
al. [sH derive a similar state from the strong coupling 
limit of a reduced two-band model. The alternative, 
extended-s wave idea gives eight nodes on the electron-FS 
sheets 11, 2^ 22 1. Our proposal is a generalization of the 
s± proposal 30l. 31 1. v4ZZ earlier proposals, however, can- 
not access possible c-axis nodal structure, which arises 
via an IB proximity effect involving dxy — d :i^2 bands 
in our work. Moreover, rigorous analysis 12 3| and the 
geometrically frustrated superexchanges 0, l2o| in FePn 
require both, Ai, A2 to be finite. This naturally leads to 
/c-dependent Aabik), as deduced experimentally pi], [s^. 
Finally, the U-SC derived here arises as an instability 
of the correlated, incoherent metal [H, I3| "normal" state 
seen in doped FePn, and naturally explains large spec- 
tral changes Q and orbital selective sharpening of the 
QP kink [l0| across Tc as manifestations of proximity to 
a Mott insulator 0, [3, 01 . 

In conclusion, we have derived the instability of the 
incoherent "normal" state of FePn Q] to an U-SC. 
Blocking of inter-band one-particle coherence in the nor- 
mal state clears the way for inter-band two-particle co- 
herence to emerge, much like in coupled LLs, giving a 
multi-orbital SC with nodes. The latter are shown to 
arise from a multiband proximity effect between the d^y 
and the c?3^2_^2 bands. Extending earlier LDA-I-DMFT 
calculations for the incoherent metal, we show how good 
quantitative agreement with the sharpening of the low- 
energy kink in the PES spectrum in the U-SC phase, as 
well as strong spectral weight transfer across the U-SC 
transition, is obtained. Further, induction of line nodes 
due to the interband proximity effect offers a reconcilia- 
tion of these attractive features with NMR and /iSR data, 
which suggest existence of line nodes in the U-SC gap. 
Taken together, these strongly support our proposal of 
an U-SC with a planar A(jfc(k) ~ (c^, -I- Cy)) + aCxCy form 
factor; i.e, a SC with co-existent, inter-site (and inter- 
orbital) s* and Sxy pair symmetries, and out-of-plane 
nodes [Ac(k) ~ Cz], at least in some of the llll-FePn 



[1] 
[2] 



[3] 
[4] 

[5] 
[6] 
[7] 

[8] 
[9] 

[10] 
[11] 



[12: 

[is: 

[14 

[is: 



[16: 

[17 

[is: 

[19 
[20 

[21 

[22: 

[23' 
[24' 
[25^ 
[26' 
[27 
[28: 

[29 
[30 

[31 
[32' 
[33' 
[34: 



Y. Kamihara et al, J. Am. Chem. Soc. 130, 3296 (2008). 

Q. Si and E. Abrahams, Phys. Rev. Lett. 101, 076401 

(2008); ibid. J. Wu et al, Phys. Rev. Lett. 101, 126401 

(2008); G. Baskaran, J. Phys. Soc. Jpn. 77, 113713 

(2008); Q. Si et al, arXiv:0901.4112 

K. Haule, et ai, Phys. Rev. Lett. 100, 226402 (2008). 

L.Craco et al, Phys. Rev. B 78, 134511 (2008); ibid M. 

S. Laad et al, Phys. Rev. B. 79, 024515 (2009). 

J. Jaroszynski et al, Phys. Rev. B 78, 174523 (2008). 

T. Imai et ai, J. Phys. Soc. Jpn. 77, 47 (2008). 

H. Luetkens et ai, Phys. Rev. L ett. 1 01, 097009 (2008); 

J.D. Fletcher et tt/.,fer Xiv:0812 .3858l 

Y.J. Uemura. [arXn^Sl 1 . 1546 ■ 

A.V. Boris et al., Phys. Rev. Lett. 102, 027001 (2009); 
ibid S.I. Mirzaei et al, arXiv:0806.2303 
L. Wray et al, Phys. Rev. B 78, 184508 (2008). 
S. Chi et al, arXiv:0812.1354 In S. Graser et al, 
arXiv:0812.0343 , and V. Mishra et al, 0901.2653, the 
pair wave function is argued to have ex-s form, and that 
disorder lifts the nodal structure. The opposite is argued 
in D. Parker et ai, Phys. Rev. B 78, 134524 (2008). 
PA. Lee et ai, Rev. Mod. Phys. 78, 17 (2006). 
P.W. Anderson, The Theory of Superconductivity in the 
High-Tc Cuprates, Princeton Univ. Press (1997). 
C. de la Cruz et al. Nature 453, 899 (2008). 
A. Nersesyan et al., in Bosonization and Strongly Cor- 
related Electronic Systems, Cambridge University Press 
(1998). 

A. Chubukov et al, Phys. Rev. B 78, 134512 (2008), 
where this issue is considered in the weak coupling limit. 

C. Xu et al, Phys. Rev. B 78, 020501 (2008); ibid C. M. 
Varma, Phys. Rev. B 73, 155113 (2006); H. Yamase and 
W. Metzner, Phys. Rev. B 75, 155117 (2007). 

M.E. Zhitomirsky and T.M. Rice, Phys. Rev. Lett. 87, 
057001 (2001). 

D. V. Khveshchenko, Phys. Rev. B 50, 380 (1994). 
M.M. Parish et al, Phys. Rev. B 78, 144514 (2008); ibid. 
T. Yildirim, Phys. Rev. Lett. 101, 057010 (2008). 

R. Fehrenbacher and M. R. Norman, Phys. Rev. Lett. 
74, 3884 (1995). 

M. Daghofer et al, Phys. Rev. Lett. 101, 237004 (2008); 

A. Moreo et aZ., arXiv:0901.3544 

Wen-Long You et al, arXiv:08 07.1493l 

A. Garg et al, Phys. Rev. B 72, 024517 (2005). 

P. Richard et al, Phys. Rev. Lett. 102, 047003 (2009). 

H. W. Ou et al, Sol. St. Commun. 148, 504 (2008). 

K. Nakamura et al, J. Phys. Soc. Jpn. 77, 093711 (2008). 
A.I. Coldea et al, Phys. Rev. Lett. 101, 216402 (2008); 
ibid. A. Carrington et a/., larXiv:0901.3976 
D. Parker et al, Phys. Rev. B 78, 134524 (2008). 

I. I. Mazin et al, Phys. Rev. Lett. 101, 057003 (2008); 
I.I. Mazin and J. Schmalian, | arXiv :0901.4790 

J. Wu and P. Phillips, arXiv:0901.3538 

L. Malone et al, arXiv:0806.3908 

Y. Ishida et al Phys. Rev. B 79 , 060503 (2009). 

J. Bauer et o?.. "a rXiv:0901.1760l 



